subroutine fme_t1fn_iefmn(fme, t1, iabij, icore, fact)
!
use mod_size
use mod_iop
use mod_ioff
use mod_orbit
implicit none
!
real*8, intent(in) :: fact
real*4, intent(in),managed :: fme(*)
!
real*4, intent(inout),managed :: t1(*), iabij(*), icore(*)
!
integer :: i, j, a, b, nsize, i0, i1, i2
integer, external :: irpdso
real*8, external :: dnrm2
!this is to do fme = fme + <ef||mn>*t(f,n)
!
i0 = 1
i1 = i0 + nvvoo*isd
i2 = i1 + nvvoo*isd
!
!===AA 
! <EF||MN>*t1(FN)+<Ef||Mn>*t1(fn)
call iexpso(iabij(ioi3(1)), icore(i0), vrta, vrta, 0, 1, popa, popa, 0, 0, 1)
call iexpso(icore(i0), icore(i1), vrta, vrta, 1, 0, popa, popa, 0, 1, 1)
call sstgenso(icore(i1), icore(i0), nsize, vrta, vrta, popa, popa, &
              icore, 1, '1324')
call VpqmnVmnrs_to_Vpqrs(icore(i0), t1(ifvo(1)), 1, 1, vrta, popa, 1, &
               vrta, popa, 1, zorb, zorb, 1, fme(ifvo(1)), 1, fact, 1.d0)
!write(6,*)'norm of fme', dnrm2(irpdso(1, vrta, popa, 1), fme(ifvo(1)), 1)
call sstgenso(iabij(ioi3(2)), icore(i0), nsize, vrta, vrtb, popa, popb, &
              icore, 1, '1324')
call VpqmnVmnrs_to_Vpqrs(icore(i0), t1(ifvo(2)), 1, 1, vrtb, popb, 1, &
               vrta, popa, 1, zorb, zorb, 1, fme(ifvo(1)), 1, fact, 1.d0)
!write(6,*)'norm of fme', dnrm2(irpdso(1, vrta, popa, 1), fme(ifvo(1)), 1)
!===BB
! <ef||mn>*t1(fn)+<eF||mN>*t1(FN)
!=<ef||mn>*t1(fn)+<Fe||Nm>*t1(FN)
call VmnpqVmnrs_to_Vpqrs(icore(i0), t1(ifvo(1)), 1, 1, vrta, popa, 1, &
               vrtb, popb, 1, zorb, zorb, 1, fme(ifvo(2)), 1, fact, 1.d0)
!
call iexpso(iabij(ioi3(3)), icore(i0), vrtb, vrtb, 0, 1, popb, popb, 0, 0, 1)
call iexpso(icore(i0), icore(i1), vrtb, vrtb, 1, 0, popb, popb, 0, 1, 1)
call sstgenso(icore(i1), icore(i0), nsize, vrtb, vrtb, popb, popb, &
              icore, 1, '1324')
call VpqmnVmnrs_to_Vpqrs(icore(i0), t1(ifvo(2)), 1, 1, vrtb, popb, 1, &
               vrtb, popb, 1, zorb, zorb, 1, fme(ifvo(2)), 1, fact, 1.d0)
!
return
end
